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Abstract 

Scattering of solitons and dark solitons by potential walls is studied 
in the nonlinear Schrodinger equation under various conditions. We in- 
vestigate the conditions under which solitons are split into two solitons 
at the potential wall. We find that a soliton can be trapped in an inter- 
space between two potential walls. A dark soliton can also be scattered 
at the potential wall. Similarly to a bright soliton, a dark soliton can pass 
through more easily the potential wall, as the width of the dark soliton 
is larger. A dark soliton can run away spontaneously from an interspace 
between the two potential walls. We also study the motion of a two- 
dimensional soliton in a two-dimensional quintic nonlinear Schrodinger 
equation. We find the coherent tunneling through a potential wall, and 
the refraction corresponding to Newton's refraction theory. 

1 Introduction 

The nonlinear Schrodinger equation is a typical soliton equation. 1 Solitons in 
optical fibers are described by the nonlinear Schrodinger equation and have been 
intensively studied for the optical communication. [21 E] Recently, solitons and 
dark solitons were observed in the Bose-Einstein condensates (BECs) . ^ [5] [7] 
Dark (bright) solitons are found in the BECs with repulsive (attractive) in- 
teraction. The nonlinear Schrodinger equation is called the Gross-Pitaevskii 
equation in the research field of BECs. In the BECs, matter waves obey the 
nonlinear Schrodinger equation with external potentials. The harmonic poten- 
tials are used to confine the BECs. Interference patterns of laser beams generate 
external periodic potentials. Solitons in external potentials have been studied 
by several authors. [51 131 ll(Jj Especially, solitons in periodic potentials (optical 
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lattices) are intensively studied. ^3 D2| However, the scattering of a soliton by 
a potential wall was not studied in detail, although it is a simple and a funda- 
mental process. The scattering problem of a quantum particle by the potential 
wall, which is described by the linear Schrodinger equation, is a typical problem 
of the quantum mechanics. In our previous letter, we studied the scattering and 
trapping of solitons by a potential wall and well. We found that the soliton be- 
haves like a classical mechanical particle under a certain condition, and evaluate 
the critical energies for the transmission. In this paper, we study a transition 
between the classical mechanical behavior and the wavy behavior more in detail. 
We will study also the scattering of a dark soliton by the potential wall. We 
will further study the scattering of a two-dimensional soliton. 

2 Scattering of bright solitons by a potential 
wall 

The nonlinear Schrodinger equation with a potential term is written as 

where 4>(x, t) is the wave function, g is the coefficient of the nonlinear term, U (x) 
is an external potential. First, we consider the case of g = — 1. It corresponds to 
the case of the attractive interaction in the problem of BECs. When U(x) = 0, 
the nonlinear Schrodinger equation with g = — 1 has a soliton solution 

<f>(x, t) = Asec\v{A{x - vt)}e tv[x - vt) - lul \ (2) 

where v — p and u> = —A 2 /2 + v 2 /2. 

In this section, we study the scattering by a simplest potential wall, which 
is expressed as 

U (x) =Uq, for xq < x < Xq + d, U (x) — 0, for x < Xq, x > Xq + d, (3) 

where d is the width of the potential wall and Uq > is the height. We study 
numerically the scattering of a soliton by the potential wall. We perform numer- 
ical simulations using the split-step Fourier method with Fourier modes 1024 
or 2048. The periodic boundary conditions are assumed, and the time step is 
0.001. Figure 1(a) displays the time evolution of a soliton passing through a 
potential wall with Uq = 0.7 and d = 1. The system size is L = 100 and the 
potential wall locates at x = xq = L/2 = 50. The initial velocity of the soliton 
is v — 1 and the amplitude is A = 1.5. The soliton is split into two pulses by 
the potential wall. The transmission coefficient is 0.488, which is calculated as 
T = J*^ , 2 1 4>\ 2 dxj \(j)\ 2 dx after the scattering. Figure 1(b) displays the trans- 
mission coefficient T as a function of the initial velocity v for three values of Uo 
for A — 1.5 and d = 1. For Uq = 0.1, the transmission coefficient changes from 
to 1 abruptly. The soliton is completely reflected below the critical value of v, 
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Figure 1: (a) Time evolution of |0(a;,t)| for Uq = 0.7,^4 = 1.5 and v = 1. 
A soliton is splitted into two pulses by the potential wall of d = 1 located at 
x = L/2 = 50. (b) Transmission coefficient T as a function of the initial velocity 
v. (c) A boundary line between the particlelike behavior and the wavy behavior. 
Under the line, a soliton behaves like a particle. The solid line is the theoretical 
estimate and the rhombi denote the critical value of Uq obtained numerically. 



and completely transmitted above the critical velocity. The soliton behaves like 
a classical mechanical particle. We showed that the critical kinetic energy de- 
pends strongly on A as E c = Uq tanh(A<i/2) using the Lagrangian formalism. |T2| 
For U — 0.7, the transmission coefficient changes smoothly from to 1. This 
behavior is characteristic of the transmission of waves or the quantum mechani- 
cal particle. Near U = 0.4, the scattering changes from the particlelike behavior 
to the wavy behavior for the amplitude A — 1.5. The transmission coefficient 
changes in a non-monotonic manner near the critical velocity, although we do 
not understand the reason yet. We discussed qualitatively the transition us- 
ing the conservation law of the energy in the previous paper. We confirm the 
conjecture with the direct numerical simulations. 

We first explain a criterion based on the conservation laws. We assume 
that a soliton is located initially at x — xo + d/2 with velocity 0. If the 
height Uq of the potential wall is sufficiently large, the soliton will be split 
into two pulses. However, the soliton cannot be split, if Uq is not so large 
owing to the attractive mutual interaction. Equation (1) has at least two in- 
variants, that is, the total mass N — J Q L \4>\ 2 dx and the total energy E = 
J Q L {1/2\4> X \ 2 - l/2\4>\ 4 + U(x)\4>\ 2 }dx. If the soliton is split into two counter- 
propagating solitons: (f>(x,t) = A'sech{A'(x - vt - L/2)}e l ^ x - vt - L / 2 ">- lu)t + 
A'sech{A'{x + vt- L/2)}e- ip ( x+vt - L / 2 '>~ luJt , the amplitude A' and the momen- 
tum p satisfy A' = A/2 and 2p 2 A' = -1/4A 3 + 2U a Ata,nh(Ad/2) owing to the 
conservation laws. If Uo < A 2 / {8 tanh(Ad/2)}, there is no real solution for p, 
that is, the splitting cannot occur and the one-pulse structure is maintained. 
Above the line Uo = A 2 / {8 tanh( Ad/ 2)}, the potential energy is used to split 
a soliton against the attractive interaction. That is, the above argument sug- 
gests that a criterion between keeping one-pulse state and splitting into two 
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Figure 2: (a) Time evolution of a soliton for Uo = 0.1, U\ = 0, A = 1 and 
v = 0.302. A soliton penetrates into a potential wall, and then is trapped 
between the double walls. The walls located at x\ — 23 and xi — 26. (b) 
Reflection coefficient R (solid line with rhombi), the trapping coefficient T p 
(dashed line with crosses) and the transmission coefficient T (dotted line with 
squares) as a function of the initial velocity v. (c) Time evolution of the peak 
position of the soliton (dashed line) and the theoretical model of eq. (6) (solid 
line). 



pulses is Uoc = A 2 /{8t&nh(Ad/2)}. To check the conjecture of the criterion, 
we have performed numerical simulations, changing the value of Uo for several 
A's, and constructed the same type plots as Fig. 1(b). We have regarded the 
critical value of U c for each A as the lowest value of Uo, at which the trans- 
mission rate T has a jump across T — 0.5. Figure 1(c) displays the critical 
values Uoc f° r A — 0.5, 1, 1.5 and 2. The dashed line is the theoretical estimate 
Uoc — A 2 / {8 tanh(^4d/2)}. The transition between the particlelike behavior to 
the wavy behavior is not so clear-cut, and so the above theoretical curve gives 
a fairly good criterion. 

The splitting of the one-soliton was caused by the perturbation at the poten- 
tial wall. It is desirable to determine the reflection and transmission coefficients 
theoretically using the inverse scattering method or other perturbative methods 
but we have not succeeded in it yet. 



3 Scattering of bright solitons by double poten- 
tial walls 

In this section, we consider the scattering of bright solitons by a potential with 
double walls: 

U(x) = Uo, for x\ < x < x\ + d, 
= —Ui, for x\ + d < x < x 2 , 
= Uq, for x 2 < x < x 2 + d, 

= 0, for x < x\, x > x 2 + d. (4) 
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Figure 3: (a) Reflection process of a soliton for Uq = 0.6, U\ = 2.5, A = 1 and 
v = 0.07. (b) Generation of a bound state at the collision for Uq = 0.6, U\ = 
2.5, A = 1 and v = 0.2. (c) Reflectivity R (solid line), trapping rate T p (dashed 
line) and transmissivity T (dotted line) as a function of v for Uq = 0.6, U\ = 2.5 
and A= 1. 

There are two potential wall at x = x\ and x 2 - Figure 2(a) displays a time 
evolution of a soliton with A = 1 and the initial velocity v = 0.302. The system 
size is 50, and the parameters of the potential are Uq = 0.1, U\ = 0, x\ = 23, 
x 2 = 26 and d = 1. The soliton passes over the first potential wall, but is 
reflected by the second wall, and then trapped in an interspace between the 
two walls. Figure 2(b) displays the reflectivity R, the transmissivity T, and 
the trapping rate T p , which are defined as R — j^ 2 5 \cf>\ 2 dx/ J Q L \(f>\ 2 dx, T p — 
Il/2-5 l^dx/ J Q L \4>\ 2 dx and T = \4>\ 2 dx/ J Q L \(f>\ 2 dx after the scattering. 

The reflectivity and the transmissivity changes from to 1 abruptly, and the 
soliton behaves like a particle. The soliton is trapped between the double walls 
in a narrow parameter region of 0.298 < v < 0.303. This trapping phenomenon 
can be described also using the Lagrangian formalism. 

The nonlinear Schrodingcr equation is rewritten with the Lagrangian form: 
d/dt{5L/6(f> t ) = 5L/5<p, where L = $ ca oQ dx{i/2{<i) t cj>* - $</>) - l/2\<p x \ 2 + 
1/2|0| 4 — U{x)\cf)\ 2 }. As is used in the previous paper, we assume a solution 
of the form 

4>{x,t) = A{t)scch{(x - ^(f))/ a ( i )} e 'P(*)(^-iW)+^(t)logco S h{(x-€(t))/a(t)}-^ ) 

(5) 

where A, a,p, £ and a change in time, but N = 2 A 2 a is conserved in time. The 
term ia(t) logcosh(x — £(t))/a(t) is a kind of chirp term. Owing to the inclusion 
of this term, the width of the soliton changes in time and breathing motion 
becomes possible. Some energy is given to the breathing motion and the kinetic 
energy may be decreased. Substitution of eq. (5) into the Lagrangian and its 
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spatial integration yields an effective Lagrangian 

L eff = N(p^-a t (l-log2)+aa t /(2a)-l/(6a 2 )-p 2 /2-a 2 /(6a 2 )+N/(6a)-U eff ), 

where U e ff = C/ /2[tanh{(a;i + d — q(t))/a} — tanh{(.xi — q{t))/a) + tanh{(x 2 + 
d — q{t))/a) — tanh{(x 2 — q(t))/a}] for U\ = 0. The equations of motion for 
p(t) , £(t) , a(t) and a(t) are 

dq dp dU e f f 

dt ^' dt dq 

da 2a da _ ( 1 + a 2 _ JV_ _ dU eff \ 

dt ~ 3a' dt V 3a 3 6a 2 da J ' ( > 

This equation is derived from a Hamiltonian H = p 2 /2 + U e ff + l/(6a 2 ) + 
a 2 /(3a 2 )+N/(6a) asdq/dt = dH/dp, dp/dt = -dH/dq, da/dt = dH/da, da/dt = 
—dH/da. The classical mechanical energy E = p 2 /2 + U e ff is not conserved 
in this equation, and the width a changes in time. If the classical mechanical 
energy E decays in time, the soliton can be trapped between the double walls. 
Figure 2(c) compares the time evolutions of q by eq. (6) with the peak position 
of the soliton by cq. (1) for the initial velocity v = p = 0.302. The trapping 
phenomenon may be qualitatively explained by eq. (6). 

When Uq and U\ become larger, a soliton is split into two or three solitons. 
We show an example of the scattering of a soliton leaving a bound state behind. 
Figure 3 displays the soliton scattering with amplitude A = 1 at Uo = 0.6 and 
U\ = 2.5. The system size is 100. Figure 3(a) displays the reflection process 
of the soliton for the initial velocity v — 0.07. For this velocity, the reflectivity 
is almost 1. But for v = 0.08, the reflectivity decreases abruptly to 0.675, 
and a bound state is created in the potential well. Figure 3(b) displays the 
reflection process leaving a bound state behind for v = 0.2. The bound state 
has a norm 0.688 and is almost stationary. The bound state is antisymmetric 
around the center, therefore, the amplitude \4>\ is at the center. The velocity 
of the reflected pulse is v ~ —0.34, and it is larger than the initial velocity. 
That is, the reflecting pulse is accelerated by leaving a bound state with low 
energy behind. Figure 3(c) displays the reflectivity R, the transmissivity T and 
the trapping rate T p as a function of the initial velocity. There is a jump of the 
reflectivity at v = 0.07. The bound state is created abruptly for v > 0.08. The 
transmissivity T increases continuously from v ~ 0.5. Near v = 1.5, a soliton is 
split into three pulses, i.e., a reflecting pulse, a bound state and the transmitted 
pulse. 

4 Scattering of dark solitons 

If the interaction is repulsive in the BECs, the nonlinear Schrodinger equation 
has a form 
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Figure 4: (a) Time evolution of \<j>(x, t)\ for U = 0.09, A = 1 and v = 3tt/50. A 
dark soliton is reflected by the potential wall of d — 1 located at x — L/2 = 25. 
(b) Time evolution of \(f>\ for Uq = 0.08,^4 = 1 and v = 37r/50. A dark soliton 
passes over the potential wall, (c) Time evolution of the minimum value of \cj>\ 
for U = 0.09 (solid line) and U = 0.08 (dashed line). 



that is, the sign of the nonlinear term becomes positive and the coefficient is as- 
sumed to be 1. A dark soliton is a special solution to this nonlinear Schrodinger 
equation. A dark soliton has a form: 

<f>(x,t) = {iB + At a nhA(x~vt)}e lk{x - vt ^ k2 ^ +kB - A2 - B ^ t 1 (8) 

where B, A and k are parameters which characterizes the dark soliton and the 
velocity v is given by k + B. For B = 0, the amplitude \<fr\ becomes at the 
center. For the amplitude does not become 0, and it is sometimes called 

a grey soliton. We consider the case of the potential given by eq. (3). The energy 
of the dark soliton is given by E — 1/2 J dx{\d4>/dx\ 2 + |0| 4 + 2U{x)\4>\ 2 } = 
-2(A 2 + B 2 )A+4/3A 3 -2kAB-k 2 A-UoA{t&nhA{xo + d-q)-taTihA(xo-q)}, 
where q is the position of the dark soliton. We investigate numerically the time 
evolution of this dark soliton. The initial value of B is set to be and the initial 
velocity is given by v = k. The width of the potential wall in cq. (3) is given by 
d = 1, and the system size 50. In actual simulations, we have investigated the 
time evolution of two dark solitons (Another dark soliton is located at q+L/2) in 
a system of double size L = 100, since we performs numerical simulations using 
the periodic boundary conditions. Owing the periodic boundary conditions, the 
initial wave number needs to be k = 2Trn/L at t = 0, where n is an integer. 
Figure 4(a) displays the reflection process of the dark soliton by the potential 
wall for U = 0.09, v = k = 37r/50, and A = 1. Figure 4(b) displays the 
transmission process for Uq — 0.08, v = 37r/50 and A = 1. The dark soliton 
approaches plane waves in the infinity, and the scattering of the plane waves by 
the potential wall always occurs, before the center of the dark soliton reaches 
the potential wall. In Fig. 4(a) and (b), the amplitude \<j>\ is slightly depressed 
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Figure 5: The critical height Uo of the potential wall, above which the dark 
soliton is reflected. The solid line is the theoretical estimate using the energy. 
The dashed line is the twice of that. 

owing to the scattering near the position x = L/2 of the potential wall, but the 
effect of the scattering of the plane waves seems to be small, when the potential 
height Uo is small. 

Figure 4(c) displays the time evolution of the minimum value of \cj>\ for 
Uq = 0.09 and 0.08. If the dark soliton keeps the form (8), the minimum value 
of \<f>\ represents \B\. The initial value of B is assumed to be 0. In the case of the 
reflection, B is decreased to -0.40, and the velocity of the reflected dark soliton is 
numerically estimated as -0.20, which is close to v = k+B ~ —0.21. The velocity 
of the reflected dark soliton is larger than the velocity v = 37r/50 ~ 0.188 of the 
initial dark soliton. Since the energy and the norm are conserved, A 2 + B 2 = 1 
and -2A + 4/3A 3 - 2kAB - k 2 A = -2/3 - k 2 . The B-value is expected as 
-0.404 from the two conservation laws if k is fixed to be the initial value, and 
it is comparable to the numerical value. When the dark soliton passes over the 
potential wall, the B- value is decreased to -0.16 near the potential wall, but it 
recovers to as the dark soliton goes away from the potential wall. 

There is a critical value of Uo for the reflection to occur. Figure 5 displays 
the critical value of Uq as a function of the initial value Aq of A above which the 
reflection is numerically observed for the initial velocity v — 37r/50. The critical 
value increases as A is decreased. It means that the dark soliton penetrates 
more easily the potential wall, as the amplitude is decreased. This effect is 
similar to the one we have already seen for bright solitons as coherent tunneling. 
The critical value may be similarly evaluated using the energy conservation law, 
if the form of the dark soliton is assumed to be given by eq. (8). At the critical 
value of Uo, the dark soliton will stop at x = x + d/2. Then, the .B-value is 
expected to be — k, since the velocity v = k + B is 0. Owing to the energy 
conservation law, -2/3 - k 2 A = -2A\A + 4/3A 3 + k 2 A - 2U Ata,-ah(Ad/2) 
is satisfied. The critical vale of Uo is evaluated as Uo = (—2A^A + A/3A 3 + 
k 2 A + k 2 A Q + 2/3Ag)/(2Atanh Ad/2). 
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Figure 6: (a) Time evolution of \<f>(x, t)\ for the double potential walls at x\ = 19 
and X2 = 30 with Uq = 0.1, U\ = 0, d = 1, A = 1 and w = 37r/50. A dark soliton 
goes out from the interspace between the double walls after several reflections, 
(b) Time evolution of the minimum point of \cf)\. (c) Time evolution of the 
minimum value of \cj>\. 



A similar result can be derived from the Lagrangian formalism. We assume 
a solution of the form 

cj){x, t) = {iB{t) + A{t) tanh A(t)(x - g (t))yMz-9«)+*(fc 2 /2+fcB-A 2 -B 2 )^ 

Substitute it to the Lagrangian 

L = J {i/2{M* - <t>* t <f>) - lmdc^/dxl 2 - 1/2|0| 4 - U(x)\cj>\ 2 }dx. 

yields an effective Lagrangian 

L eff = -2BAq t - 2A(kq t - uj) + 2(A 2 + B 2 )A - 4/3A 3 + Ak 2 + 2ABk 

+U a A{ta,-nh{A(xQ + d - q)} - tanh{A(x - q)}, 

where ui = k 2 /2 + kB and A 2 + B 2 = A\ = const. We can obtain equations of 
motion from the effective Lagrangian as 

d(BA + Ak) = A dU e/f 
dt dq 

d JL k + B+ Z d(AU eJf ) 

dt " + + B^-A^ + kB OA ' W 

where U e fj = l/2U {ta,nh{A(x + d — q)} — tanh{A(x — q)}]. Similar critical 
value of Uq can be evaluated using this equation, which is almost the same as the 
one evaluated using the energy. However, the critical value of Uq estimated by 
the above theory is close to half of the critical value by numerical simulations as 
shown in Fig. 5. The above equations of motion are approximated to d 2 q/dt 2 = 
—dU e ff/dq when v and Uq are sufficiently small and then \B\ is small, which 
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Figure 7: (a) Amplitude profile f(r) = \<f>(x, y, t)\ for the two-dimensional 
quintic nonlinear Schrodinger equation with c = 0.1. (b) Time evolution of 
\4>(x, L/2)\ for the initial velocity k x — 0.19. (c) Time evolution of \<p(x, L/2)\ 
for the initial velocity k x = 0.24. 

is the equations of motion for a classical particle of mass 1. On the other 
hand, Busch and Anglin derived equations of motion of the form d 2 q/dt 2 = 
— (l/2)dU(q)/dq when the potential U(x) changes very slowly compared to the 
soliton scale. [13] using a multiple time scale boundary layer theory. In their 
theory, the background amplitude A 2 + B 2 and the velocity are assumed to 
change slowly. The resultant effective mass of the dark soliton becomes 2. 
Considering the factor 1/2, it is reasonable that the numerical critical value of 
Uq is twice of the theoretical estimate. That is, our evaluation based on the 
energy and the Lagrangian fails to estimate the critical value of Uq. Probably, 
it is necessary to involve the effect of the variation of the background component 
A, B, and k to the Lagrangian formalism, although we have not succeeded yet. 

For a bright soliton, the kinetic energy is somewhat decreased and can be 
trapped in an interspace between double potential walls, even if the soliton is 
initially outside of the potential walls. On the other hand, the dark soliton 
can get some kinetic energy at the scattering and run away from the interspace 
between the double walls, because the dark soliton with smaller B-value is 
energetically disadvantageous. To show the detrapping, we have used the double 
potential walls of eq. (4). The parameters are Uq — 0.1, U\ = 0, Xi = 19 and 
x% = 30 and d = 1. The initial amplitude A is 1 and the initial velocity is 
v = 37r/50, and the initial position is q — 25. Figure 6(a) displays the time 
evolution of the dark soliton. Figure 6(b) displays the time evolution of the 
minimum point of \<f>\. Figure 6(c) displays the time evolution of the minimum 
value of \cf>\. The maximum of the |-B|-value increases slowly each time the dark 
soliton is reflected by the second wall, and finally the dark soliton overcomes 
the potential wall. This detrapping phenomenon can be described by eq. (9), 
although it takes more time to overcome the potential walls in eq. (9). We need 
to improve the equations of motion to describe the detrapping phenomenon 
more quantitatively. 
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Figure 8: (a) Time evolutions of X p and Y p (t), where (X p ,Y p ) is the peak 
position of the two-dimensional soliton. (b) Trajectory (X p (t),Y p (t)) of the 
peak position in the x — y plane. The initial velocity is (0.3,0.1). The dashed 
line denotes a boundary x — L/2 of the stepwise potential. 

5 Scattering and refraction of two-dimensional 
solitons 

The standard two dimensional nonlinear Schrodinger equation does not have 
soliton type solutions. Even if a localized initial condition is assumed, it may 
extend and changes into radiations or leads to collapse, depending on the initial 
conditions. However, two dimensional solitonlikc solutions exist in the quin- 
tic nonlinear Schrodinger equation. [3] The two-dimensional quintic nonlinear 
Schrodinger equation is not an integrable system, so it may not be suitable to 
call the localized solution a soliton, but we say it a two-dimensional soliton in 
this paper by analogy with the one-dimensional system. We study the scatter- 
ing and the refraction of the two-dimensional soliton using the quintic nonlinear 
Schrodinger equation: 

* W = ~ + + U (X ' V)(t> - (10) 

Equation (10) has an axisymmetric localized solution for U(x, y) — 0. The 
localized solution has a form <j>(x,y) = f(r)e~ ltJ,t , where r is the distance from 
the center of the axisymmetric localized solution. The amplitude /(r) needs to 
satisfy 

f" + r-\f> = -2{n + f 2 -cf}f. (11) 

Figure 7(a) is a numerically obtained solution to the above nonlinear eigenvalue 
problem for /(0) = 2 and c = 0.1. The amplitude /(0) of the two-dimensional 
soliton is a free parameter, as is the one-dimensional soliton. We have checked 
that the axisymmetric solution is stable by direct numerical simulations of 
eq. (10). We can construct a propagating soliton with velocity (v x ,v y ) = (k x , k y ) 
using the solution of eq. (11) as <j> — f(x — v x t, y ~ Vy t^ e ik ^ x +' lk yy-nt *_ We have 
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performed numerical simulations of the quintic nonlinear Schrodinger equation 
with a potential wall of U(x,y) = U = 0.04 for L/2 - 0.5 < x < L/2 + 0.5. 
That is, the potential wall is parallel to the y-axis and has a width d = 1. 
The system size is 25 x 25, and the initial wavenumber is (k x ,0), that is, the 
initial velocity is perpendicular to the potential wall. Figure 7(b) and (c) dis- 
play the time evolutions of \<fi(x,L/2)\ for k x — 0.19 and 0.24 at y = L/2. 
We have checked that the form of the two-dimensionally localized structure 
is maintained at the scattering. The two-dimensional soliton is reflected for 
k x = 0.19 and it has passed through the potential wall for k y = 0.24. The 
critical value for the transition is k x = 0.215 in our numerical simulation. 
The Lagrangian for the quintic nonlinear Schrodinger equation is given by 
L = dxdy{i/2(M* ~ <f>* t 4>) - 1/2|V0| 2 + 1/2|0| 4 - c/3|0| 6 - U(x, y)\<j>\ 2 }. 
If the solution is assumed to be <p(x, y) = f(x — q{t),y)e lp< ^ x ~ q ^~ %tlt , the equa- 
tion of motion for the center of mass is written using the Lagrangian formalism 
d/dt(dL e ff/dq) = dL e ff/dq, d/dt(dL e ff /dp) = 8L e ff/dp as 



d*q = dUeff 

dt 2 dq ' 



(12) 



where U e ff — J J f 2 (x — q,y)U(x,y)dxdy/ J 2-rrr f (r) 2 dr . For the potential wall 
with height Uq and the width d, the critical velocity v c is estimated from the 
relation 

1 2 _ Iq /2 ^rf 2 dr + ^ 2 r4sm- 1 {d/(2r)}f 2 dr 

2 Vc ~ J °°2nrf 2 dr ' ( 3 > 

For d = 1 and Uo — 0.04, the critical initial velocity is 0.214 for f(r) shown in 
Fig. 7(a), which is consistent with the direct numerical simulation. The critical 
energy is sufficiently smaller than the potential height Uo similarly to the one- 
dimensional case. It might also be interpreted as the coherent tunneling. 

We have performed a numerical simulation for the steplike potential: U(x) = 
for x < L/2 and U(x) — Uo for x > L/2. The initial velocity vector is 
(k x , k y ) = (0.3, 0.1), and the potential height is Uo = 0.04. Figure 8(a) displays 
time evolutions of the peak position (X p ,Y p ) of the two-dimensional soliton. 
The velocity in the x-direction changes from 0.3 to 0.1, and the velocity in the 
y-direction does not change. The velocity v x — 0.1 in the x-direction after the 
scattering satisfies the energy conservation law l/2v 2 + Uo = 1/2 • 0.3 2 . As a 
result, the trajectory of the two-dimensional soliton is refracted, which is shown 
in Fig. 8(b). The angle of the refraction is nearly 7r/4, and it is larger than the 
angle of the incidence, although the velocity becomes smaller after the refraction 
for x > L/2. This phenomenon corresponds to Newton's particle theory for the 
refraction. The soliton behaves like a particle, and therefore Newton's theory is 
revived. 
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6 Summary 



We have studied the scattering of bright solitons and dark solitons by the 
potential walls. We have confirmed by direct numerical simulations that the 
boundary line of the particlelike behavior and the wavy behavior is given by 
Uqc = A 2 / {8 tanh(^4<i/2)}, which was expected in our previous paper. We 
have found that a bright soliton can be trapped between the double potentials. 
When Uq is small, the soliton behaves like a particle and the soliton is com- 
pletely trapped in a narrow parameter region between the reflection and the 
transmission. When Uq is large and U\ is large, there appears a bound state 
in the interspace between the double walls. At the collision of a soliton to the 
left wall, a part of the soliton is changed into the bound state and the rest 
is reflected. The velocity of the reflected pulse can be faster than the initial 
velocity. 

We have also studied the scattering of dark solitons. When the dark soliton is 
reflected, the B-value is decreased from to a certain negative value, that is, the 
dark soliton is changed into a grey soliton. We have evaluated the critical value 
of Uo between the reflection and the transmission. A simple theoretical estimate 
is about half of the numerical values. It may be because the background density 
and velocity vary in space by the potential wall. The calculation by Busch and 
Anglin is based on the approximation of the slowly varying potential. For the 
potential walls whose width are comparable to or smaller than the soliton's scale, 
their approximation does not seem to be applicable directly, but the factor of 
1/2 may be related to the effect. We need to improve the theoretical estimate 
by considering the effect. We have also found that a dark soliton cannot be 
trapped in an interspace between the double walls, and runs naturally away 
from the trapping region. 

We have further studied the scattering and the refraction of a two-dimensional 
soliton in the quintic nonlinear Schrodinger equation. We have found that the 
coherent tunneling occurs also for the two dimensional soliton. We have found 
a phenomenon corresponding to Newton's theory for the refraction. 
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